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ABSTRACT 

We study effects of external ionization by cosmic rays (CRs) and X-rays on the ther- 
mal evolution of primordial clouds under strong far-ultraviolet (FUV) radiation. A 
strong FUV radiation field photodissociates H 2 and quenches its cooling. Even in such 
an environment, a massive cloud with the virial temperature > 10 4 K can contract 
isothermally at 8000 K by the hydrogen Lyman a cooling. This cloud collapses mono- 
lithically without fragmentation, and a supermassive star (> 10 5 Mq) is believed to 
form at the center, which eventually evolves to a supermassive black hole (SMBH). 
However, candidates of FUV sources, including star-forming galaxies, are probably 
sources of strong CRs and X-rays, as well. We find that the external ionization pro- 
motes H2 production and elevates the threshold FUV intensity J CT n needed for the 
SMBH formation for CR energy density Ucr > f0~ 14 erg cm -3 or X-ray intensity 
Jx > 10~ 24 erg s _1 cm~ 2 sr _1 Hz" 1 at 1 keV. The critical FUV flux increases as 

J cr it oc Ucr 1 ^ 2 (oc Jx 1 ^ 2 ) in the high CR (X-ray, respectively) limit. With the same 
value of FUV intensity at the Lyman limit (13.6 eV), the H~ photodissociation rate, 
with threshold of 0.755 eV, increases and the H2 abundance decreases with decreasing 
effective temperature of the FUV sources T* . The lower value of T* thus results in the 
lower critical FUV flux </ cr it at the Lyman limit. Using an empirical relation between 
intensities of FUV and CRs/X-rays from nearby star- forming galaxies, we find that 
external ionization effect remarkably enhances the critical FUV flux for sources with 
T* as high as 10 5 K and composed of stars with < 100 Mq to a level that is not 
realized in any halo. This indicates that, to induce SMBH formation, the FUV sources 
must be either Pop II /I galaxies with low brightness temperature (T* ~ 10 4 K), Pop 
III galaxies (T* ~ 10 5 K) with a very top-heavy IMF, or Pop III galaxies too young to 
harbor sources of CRs/X-rays, for example, supernova remnants or high- mass X-ray 
binaries. 
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1 INTRODUCTION 

The origin of supermassive black holes (SMBHs), ubiqui- 
tously existing at the centers of present-day galaxies, still 
remains a mystery. Among them, the existence of the lumi- 
nous SDSS quasars with estimated SMBH masses of a few 
times 10 9 Mq by their luminosities > 10 47 erg s _1 at red- 
shift z > 6, when the age of Universe is less than 1 Gyr, 
(e.g., Fan 2006) poses serious constraints on their formation 
and evolution scenarios. Several authors have studied mod- 
els for the SMBH growth by merger and gas accretion (e.g., 
Haiman & Loeb 200f ; Volonteri, Haardt & Madau 2003; Li 
et al. 2007), starting from stellar-mass black holes (BHs), 
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which are remnants of first generation of stars. However, 
various negative feedbacks prohibit their rapid growth. 

In a merger event, the BH experiences strong recoil, 
with typical velocity > fOO km s _1 , owing to gravitational 
wave emission (Herrmann et al. 2007; Koppitz et al. 2007). 
The recoil velocity depends on the spin configuration of the 
merging pair and can reach as high as ~ 4000 km s" 1 for the 
anti-parallel spins (Campanelli et al. 2007). In high-redshift 
universe, typical halos tend to be low-mass and their escape 
velocity (< fO km s _1 ) is smaller than the recoil velocity. If 
BHs are ejected from their host halos in merging events, the 
BH growth process must repeat again from scratch. On the 
other hand, in the gas-accretion scenario, the accretion rate 
onto the BH is at most similar to the Eddington limit. In this 
case, it takes ~ f Gyr to form a SMBH with 10 9 Mq from 
a BH with fOO Mq. Although this is marginally consistent 
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with the existence of SMBHs at z > 6, this requires constant 
supply of gas at the Eddington rate. However, recent hydro- 
dynamical simulations demonstrate that radiative feedback 
from the the accreting black hole prevents such efficient ac- 
cretion (Milosavljevic, Couch & Bromm 2009; Alvarez, Wise 
& Abel 2009). 

We here explore another possibility that the seed BHs 
form via gravitational collapse of supermassive stars (SMSs) 
and thus are as massive as > 10 s Mq from the beginning. 
Evolution and general relativistic instability of the super- 
massive stars have been investigated by many authors (e.g., 
Chandrasekhar 1964a, b; Zeldovich & Novikov 1971; Shapiro 
& Teukolsky 1983). Among them, Shibata & Shapiro (2002) 
calculated the collapse of a rotating supermassive star into 
a SMBH and found that most of the mass originally in the 
SMS is eventually locked in the BH (M BH ^ 0.9 M S ms). 

The biggest challenge in this scenario is how such super- 
massive stars are formed. A promising pathway is through 
the collapse of a cloud in a halo under an extremely strong 
far-ultraviolet (FUV) radiation field. The FUV radiation in 
the Lyman and Werner bands (11.2—13.6 eV) dissociates Eb 
molecules, which are important coolants in a primordial gas 
in the low-temperature (< 8000 K) regime. Under a FUV 
radiation field exceeding a critical value, J cr it ^ 300 (10 ) 
in the unit of 10 -21 erg s _1 cm' 2 sr _1 Hz -1 , for the stellar- 
type radiation with brightness temperature T» = 10 4 K (10 5 
K, respectively), the H2 formation is suppressed all the way 
during the protostellar collapse and then its cooling is totally 
quenched (Omukai 2001; Shang, Bryan & Haiman 2010). 
In this case, the collapse proceeds almost isothermally at 
~ 8000 K by hydrogen Lyman a and two-photon emission, 
and HP bound-free emission. Unlike the case with weaker 
radiation field, where the H2 line cooling induces fragmen- 
tation to clumps of ~ 10 3 Mq, such a cloud eludes frag- 
mentation and collapses directly to a supermassive star of 
> 10 5 Mq (Bromm & Loeb 2003; Regan & Haehnelt 2009a, 
b), and collapses eventually into the seed BH for the subse- 
quent growth to a SMBH. 

The critical FUV flux for the suppression of the H2 
cooling is far above the expected mean value of the back- 
ground radiation. Studying the spatial correlation between 
halos, Dijkstra et al. (2008) derived the probability distribu- 
tion of the intensity of Lyman- Werner flux J^i* ( m units of 
10~ 21 erg s _1 cm~ 2 sr _1 Hz -1 ) incident on halos with mass 
~ 10 Mq collapsing at redshift z ~ 10. The estimated mean 
Lyman- Werner background is = 40, far below the crit- 
ical value. This means that such intense FUV fields are only 
realized in neighborhoods of strong FUV sources. They also 
estimated the fraction of halos bathed in radiation fields ex- 
ceeding the threshold value Jcrit ■ This fraction is 10 for 
Jcrit = 10 3 and decreases exponentially with increasing J cr i t 
(> 2 x 10 3 ). A small difference in J cr i t significantly affects 
the number of halos bearing supermassive stars, and thus 
correct knowledge of the critical FUV flux is crucial in esti- 
mating the number of seed BHs. 

So far, the critical FUV flux J cr i t has been studied only 
in cases where the incident radiation consists of components 
below the Lyman limit (Omukai 2001; Bromm & Leob 2003; 
Shang et al. 2010). However, strong FUV sources, i.e., ac- 
tively star-forming galaxies, are expected to contain a large 
number of massive stars and possibly some mini-quasars. 
Since massive stars end their lives as supernovae and leave 



the remnants, where cosmic rays (CRs) are accelerated and 
X-ray photons are produced, sources of strong FUV radia- 
tion can also be those of CRs and X-ray photons. Similarly, 
mini-quasars and high-mass X-ray binaries emit soft X-ray 
radiation. The incident flux is thus expected to have such 
high-energy components. If present, CRs and X-rays en- 
hance the ionization degree of gas and increase the amount 
of H2, which is formed by electron-catalysed reactions. In 
fact, Haiman, Rees, & Loeb (1996) and Glover & Brand 
(2003), who studied the condition for virialized minihalos 
under both FUV and X-ray irradiation to be able to cool 
via H2 line emission, found that the FUV photodissociation 
effect is somewhat alleviated by the X-ray ionization. 

In this paper, we consider similar effects of X-ray and 
CR ionization on more massive halos with the virial tem- 
perature > 10 4 K. Even without H2 cooling, the clouds in 
such halos start collapsing via H atomic cooling. If the FUV 
field is below some critical value Jcrit, however, abundant 
H2 molecules form and fragmentation of the clouds occurs 
at some high density during the collapse by the H2 cooling. 
With the enhanced H2 fraction by X-ray and CR ioniza- 
tion, the more FUV flux is necessary to quench its cooling, 
thereby boosting the critical value J cr it- This results in fewer 
number of seed BHs since the isothermal collapse via atomic 
cooling is required to form SMSs. Our aim is to see the de- 
pendence of the critical FUV flux J cr it on the amount of CRs 
as well as X-rays. By this way, we try to discuss the nature 
of the FUV sources needed to induce SMS formation. 

The organization of this paper is as follows. In Section 
2, we explain our model of calculation, including chemical 
reactions and effects of FUV, CRs/X-rays. We present our 
results for the evolution of the clouds irradiated by FUV 
radiation and ionized by CRs/X-rays in Section 3, and give 
analytic interpretation for the obtained critical FUV flux 
in Section 4. In Section 5, we use an empirical relation be- 
tween intensities of FUV and CRs/X-rays from nearby star- 
forming galaxies and speculate conditions to induce the su- 
permassive star formation. Finally, we summarize this paper 
and give some discussions in Section 6. 



2 MODEL 

2.1 thermal evolution 

We consider thermal evolution of a metal-free cloud in a 
moderately massive halo with virial temperature > 10 4 K 
irradiated by a FUV field and simultaneously by either CRs 
or X-rays. We suppose the radiation sources are extragalac- 
tic background or nearby other halos, rather than local ones 
in the same halo. In star-burst galaxies, although the inter- 
nal radiation field can be very high, the interstellar media 
are expected to be metal-enriched. Since even a low level of 
metal enrichment (> 10 -5 Zq) induces fragmentation and 
prevents the SMS formation (Omukai, Schneider & Haiman 
2008), we do not consider here the cases of internal sources. 

The evolution during the gravitational collapse is calcu- 
lated by a one-zone model (e.g., Omukai 2001). The actual 
hydrodynamical collapse of self-gravitating clouds is well de- 
scribed by the Penston-Larson self-similar solution (Penston 
1969; Larson 1969). In this solution, the clouds consist of 
two spatial parts: the flat central core, whose size is about 
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the Jeans length Aj and where number density of hydro- 
gen nuclei tih, temperature T, chemical concentrations of 
species i, y(i), etc. are nearly homogeneous, and the enve- 
lope, where the density profile obeys the power law with ra- 
dius as p oc r~ 2 . The central density increases approximately 
at the free-fall rate. In our one-zone model, all the physical 
variables such as i%h, T, etc. are intended to indicate those at 
the center of the core. It has been confirmed that the thermal 
evolution of primordial clouds under strong FUV fields but 
without X-ray/CR ionization by the one-zone model repro- 
duces very well that by the full three-dimensional calculation 
(Shang et al. 2010). 

Our main interest in this paper is put on the cases with 
extremely strong FUV radiation, where H2 is totally dis- 
sociated and its cooling is quenched. To start gravitational 
collapse, such halos must be as massive as with its virial tem- 
perature > 10 4 K, where the H atomic cooling is effective. 
This requires the total mass of the halo being > 2 x 10 7 Mq, 
corresponding to > 4 x 10 Mq of gas for the virialization 
epoch 2 v ir — 10. Since the pressure effect is not important 
in such massive halos, the baryonic density in the clouds p 
is assumed to increase with the free-fall timescale: 

dp = p_ 

dt tff ' 



(1) 



where tg = y /, 3n/32Gp. In evaluating ts, we neglect the 
contribution from the dark matter. The dark-matter grav- 
ity dominates that by the gas and shortens the collapse 
timescale only just after the virialization, where the tem- 
perature increases adiabatically in any case. We confirmed 
from experiments that the omission of dark-matter grav- 
ity does not cause any noticeable deviations in temperature 
evolution as a function of density. 

The temperature evolution is calculated by solving the 
energy equation: 

de d /1\ A nct 

di = - p dA-p)-— (2) 

where e is the internal energy per unit mass, 

: An7 " 0) 



Tad - 1 /WnH ' 



p is the pressure, 



P - 



pk B T 
prriH ' 



(4) 



p is the mean molecular weight, fee is the Boltzmann con- 
stant, 7 a< j is the ratio of specific heat and A nct is the net 
cooling rate per unit volume. We set p — 1.2 and 7 a( j = 5/3 
in our calculation. The net cooling rate consists of the fol- 
lowing contributions: 



Ah + Ah 2 + Ahd + A chom — Tor — T> 



(5) 



where Ah, Ah 2 and Ahd are the radiative cooling rates by H 
atom, H2 molecule and HD molecule, respectively, A c h cm is 
the net cooling rate associated with chemical reactions, and 
Tor and Tx are the heating rates by CRs and X-rays, re- 
spectively. In the metal-free gas, H cooling by Ly a and two- 
photon emissions is effective at temperature > 8000 K, while 
H2 molecule is an important coolant at lower temperature. 
At even lower temperature < 150 K, HD molecule becomes 
the dominant coolant. Such low temperature environment 
is realized for example in the case of ionization by intense 



CRs or X-rays, which promotes H2 formation via electron 
catalysed reactions, and thereby lowering the temperature 
below ~ 150 K. The H atomic cooling rate is calculated by 
solving the level populations as in Omukai (2001), while, for 
the H2 and HD cooling rate, we adopt the fitting formulae 
by Galli & Palla (1998). In our calculation, we neglect the 
radiative cooling by He species, which is only effective for 
T > 10 5 K. The methods of calculation for heating rates 
by CRs and X-rays are described in Section 2.4.2 and 2.4.3, 
respectively. 

When the density increases and the cloud becomes op- 
tically thick, the intensities of external radiation, i.e., FUV, 
CRs and X-rays, are attenuated before reaching the center 
and also the radiative cooling by e.g. H Ly a and H2 becomes 
ineffective by the self-absorption. Wc assume the radius of 
the central core R c is a half of the Jeans length Aj; 



Rc - T - 2 



TTk B T 

GpprriH ' 



(6) 



Since we focus on the center of collapsing clouds, the column 
number density of i-th chemical species is 



N(i) = y(i)nnRc- 



(7) 



We use this column density in considering the optical-depth 
effect on the incident radiation intensities, as well as on the 
radiative cooling rates. 

2.2 chemistry 

We consider the chemical reactions in primordial gas among 
the following 14 species: H, H 2 , e", H + , H 4 ", H", D, HD, 
D+, HD+, D", He, He+, and He ++ . The included reactions 
are summarized in Table 1. The following points are worth 
noting. We adopt the H2 collisional dissociation rate (reac- 
tion 7) of Martin, Schwarz & Mandy (1996), which is ten 
times larger than the older rate of Shapiro & Kang (1987) 
used by Omukai (2001) around ~ 10 3 cm~ 3 . Shang et al. 
(2010) found that this difference affects the value of the crit- 
ical FUV flux Jcrit by about an order of magnitude. For the 
photoionization of H (reactions 21 — 25) and the photodis- 
sociation of H^ (reaction 20), those from the excited levels 
are also included as in Omukai (2001). 

2.3 Initial condition 

The calculation starts at nu — 0.1 cm~ 3 , with initial tem- 
perature T — 160 K, appropriate for halos virializing at 
Zvir = 10, whose turnaround epoch corresponds to 2turn — 17 
(Omukai et al. 2008). The early stages at densities tih = 0.1 
cm~ 3 and 1.0 cm~ 3 correspond to redshifts of z ~ 15.5 
and 11.5, respectively. The initial concentrations of electron, 
y(e~) = 2 x 10 -4 , and H 2 molecule, y(U 2 ) = 2 x 10~ 6 , arc 
taken from the values at ~ 0.1 cm~ 3 of Bromm & Loeb 
(2003). We set the initial He concentration at y(He) = 0.08, 
which corresponds to the He mass fraction, Yho = 0.24. 



2.4 Incident FUV radiation, CRs, and X-rays 

2.4.1 FUV radiation 

The incident FUV radiation field Jfuv(^) is assumed to 
have a diluted thermal spectrum, i.e., Jfuv(^) oc _B„(T»), 
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Table 1. Chemical Reactions 



Number 


Reaction 


Reference 


H collisional reactions 


1 


H + e- -> H+ + 2e" 


4 


2 


H+ + e~ ->■ H + 7 


3* 


3 


H + e~ -> H~ + 7 


2 


1 


H- + H -> H 2 + c 


2 


5 


H + H+ -> H+ + 7 


2 


6 


H^4- H -> H2 + H+ 


2 


7 


H 2 + H -> 3H 


4 


8 


H 2 + H+ -v H+ + H 


2 


9 


H 2 4- e- -> 2H + e~ 


2 


10 


H~ + e" -)H + 2e" 


1 


11 


H- 4- H+ -> 2H 


2 


12 


H~ + H+ — ► 4- e~ 


2 


13 


4- e- -> 2H 


2 


14 


H+ 4- H- -> H 2 4- H 


1 


15 


3H -J- H 2 4- H 


5 


16 


2H 4- H 2 -> 2H 2 


5 


17 


2H 2 -> 2H 4- H 2 


5 


photo-, and CR reactions 


18 


H 2 4- 7 ->■ 2H 


6 


19 


H"47->H4e" 


7 


20 


H 2 + 47->H4H+ 


8 


21 - 25 


H(n) + 7 -> e - + H+ (n = 1 - 5) 


9 


26 


H 4- CR -* e" 4- H+ 


10 


D reactions 


27 


D+ + e" -) D + 7 


11 


28 


D4H+^D+4H 


12 


29 


D+ 4- H -> H+ 4- D 


12 


30 


D 4- H -> HD 4- 7 


11 


31 


D 4 H2 -> H + HD 


11 


32 


HD+ 4- H -s. H+ 4- HD 


11 


33 


D+ + H 2 -> H+ 4- HD 


13 


34 


HD + H -> H 2 4- D 


11 


35 


HD 4- H+ -> H 2 4- D+ 


13 


36 


D + H+-> HD+ 4- 7 


11 


37 


D+ 4- H -> HD+ 4- 7 


11 


38 


HD+ 4e'->H + D 


11 


39 


D + e"^D"+7 


11 


40 


D+ 4- D~ -> 2D 


11 


41 


H+ 4- D~ -» D + H 


11 


42 


H~ 4- D — > H 4- D~ 


11 


43 


D"4H-)D + H- 


11 


44 


D~ 4- H — > HD 4- e~ 


11 


45 


HD 4- 7 -+ H 4- D 


6 


He reactions 


46 


He 4- e~ -> He+ 4- 2e" 


1 


47 


He+ 4- e~ ->■ He 4- 7 


1 


48 


He+ + e - -* He++ 4- 2e _ 


1 


49 


Hc++ + e~ -> He+ 4- H+ 4-7 


1 


50 


He+ 4- H -> He 4- H+ 4- 7 


14 


51 


He 4- H+ -> Hc+ 4- H 


14 


52 


He 4- 7 ->• He+ 4- e - 


1 



(1) Abel et al. (1997); (2) Galli & Palla (1998); (3) Ferland et 
al. (1992; Case B); (4) Martin et al. (1996); (5) Palla, Salpeter 
& Stahler (1983); (6) Wolcott-Green & Haiman (2011); (7) John 
(1988); (8) Stancil (1994); (9) Rybicki & Lightman (1979); (10) 
Stacy & Bromm (2007); (11) Nakamura & Umemura (2002); (12) 
Savin (2002); (13) Galli & Palla (2002); (14) Yoshida et al. (2006). 




column density Nh (cm- 2 ) 



Figure 1. The cosmic-ray ionization rate Ccr f° r ^CR = 10~ 15 
erg cm -3 and e m i n = 10 6 eV, and e ma x = 10 15 eV as a function 
of H column density TVh- In the high column-density limit, the 
rate is approximated as fcR oc N^ ' 75 . 



with brightness temperature T* = 10 4 or 10 5 K, correspond- 
ing to that from assembly of metal-enriched stars or mas- 
sive Pop III stars, respectively. The normalization of the 
intensity is set at the Lyman-limit frequency — 13.6 eV 
by J21 = JfuvM/KT 21 erg s~ x cm' 2 sr _1 Hz -1 . For the 
same value of J21 at the Lyman limit, the rate of H 2 pho- 
todissociation, which proceeds via absorption of photons of 
11.2 — 13.6 eV does not vary so much for different values 
of T„. On the other hand, H~ photodissociation rate, whose 
threshold energy (0.755 eV) is far below the Lyman limit, 
is significantly affected with change of T* even with the 
same J 2 i. For example, this rate is 2 x 10 4 times higher 
for T» = 10 4 K than that for T* = 10 5 K with the same J 2i . 
Since H~ is the intermediary in the H 2 -forming reaction (re- 
actions 3 and 4), the H 2 concentration depends sensitively 
on T* (see Section 4). 

The H 2 and HD are self-shielded against photodissocia- 
tion for their column densities > 10 13 cm -2 . The HD shield- 
ing is also contributed by H and H 2 . We use the shielding 
factors by Wolcott-Green & Haiman (2011). 



2.4-2 Cosmic Rays 

For the incident CR flux, we follow the treatment by 
Stacy & Bromm (2007). The CR energy distribution 
dncYi/de (cm -3 eV -1 ) is assumed to obey a power-law spec- 
trum with index —2, the value by the diffusive shock acceler- 
ation by a strong shock (e.g. Bell 1978) from the minimum 
energy e m i n = 10 6 eV to the maximum e max = 10 15 eV. 
After incident to the gas, CRs propagate inward with los- 
ing their energies by ionization up to the penetration depth 
Dp at a rate (de/dt) ion, for which we use the expression by 
Schlickeiser (2002). The heating rate by CRs at depth D is 



r C n(D) = 




where Ae = 50 eV is the approximate kinetic energy of a 
CR particle lost upon each scattering (Spitzer & Tomasko 
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1968), and the ionization rate is 



Ccr(O) = 



r C R(D) 



(9) 



where £a 1C at — 6 eV is the energy deposited as heat per 
ionization in a neutral medium (Spitzer & Scott 1969; Shull 
& van Steenberg 1985). Note that the CR ionization rate is 
related with the total CR energy density, 



dncR , 



by 



Ccr = 1.7 x 10 



-18 rr -1 

Uis s , 



(10) 



(11) 



in the low column density limit, where Uis = Ucn/W~ 15 
erg cm -3 . In Fig. [TJ we also plot the CR ionization rate 
Ccr as a function of the column density Nh- This can be 
approximated as 



6.8 x 10" 



-3/4 

T 



10 3 cm- 3 8000 K 



-3/8 



\ (12) 



for column density higher than > 10 22 cm -2 . 

The CR ionization rate in the Milky Way has been 
estimated by many authors, including Ccr. ~ 4 x 10 -16 
s _1 (Hayakawa, Nishimura & Takayanagi 1961), 6.8 x 
10" 18 s _1 < Ccr < 1-2 x 10~ 15 s" 1 (Spitzer & Tomasko 
1968) and Ccr ^ 3 x 10 -17 s" 1 (Webber 1998). By recent 
observation of HjJ" in the interstellar medium, the average 
H2 ionization by CRs is evaluated Ccr,h 2 — 4 x 10 -16 s _1 
(McCall et al. 2003; Indriolo et al. 2007), which is translated 
to a rather high H ionization rate of Ccr — 2.6 x 10 -16 s _1 . 
The CR intensity at high redshift is totally uncertain. The 
CR energy density is set by a balance between the injec- 
tion and the diffusive leakage. On one hand, star forma- 
tion in young galaxies can be more active than in the Milky 
Way. The CR injection rate being proportional to the star 
formation rate, the gas in the neighborhood of such galax- 
ies is subject to intense CRs. On the other hand, magnetic 
fields are expected to be much weaker in young galaxies, 
which results in the longer Larmor radii of CRs and eas- 
ier leakage from the galaxies. This may result in weaker 
CR density if the CR sources are in the same halo. Con- 
sidering those uncertainties, we here regard the CR en- 
ergy density as a free parameter and calculate the cases for 
10~ 3 ^ Uis ^ 10 4 , which corresponds to the CR ionization 
rate of 10 -21 s _1 < Ccr i$ 10~ 14 s -1 in the low density case. 



2.4.3 X-rays 

Since the ionization cross sections of hydrogen and helium 
fall as an(u) oc v~ 3 and <TH e (y) oc v~ 2 , respectively, towards 
higher energy photons, the soft X-ray (2 — 10 keV) photons 
reach far longer distance from the sources than ionizing UV 
photons. Many sources contribute to X-rays at a given point, 
i.e., an extragalactic X-ray background is built up (Haiman, 
Rees & Loeb 1997). 

Following Glover & Brand (2003), we assume the inci- 
dent X-ray background having a power-law spectrum with 



index —1.5, 

JxM — Jx 21 x 10~ 21 ( — ') erg s _1 cm -2 sr _1 Hz -1 , 
Wo/ 

(13) 

where hvo — 1 keV. We consider the X-ray ionization of 
both H and He atoms. The heating rate by X-rays is 



Tx = Tx,H + Tx.He 



(14) 



where 



!\ = / ^^ e^^iW^ (i = H, He), (15) 



hv 



the optical depth is given by 



Tv = NhCTh(v) + A^HeCTHe(^), 



(16) 



Eh,i is the energy deposited as heat for each ionization pro- 
cess, and given by a formula by Wolfire et al. (1995). The 
ionization rate is also expressed as 



Cx 



X,p 



4ttJ x (^) 
hv 



)dv, 



(17) 



where the subscript p represents primary ionization by X- 
rays. The primary electron's energy is deposited not only in 
heating but also in the secondary ionization. Since the en- 
ergy of the primary electron is far larger than the ionization 
threshold, the secondary ionization is more effective than 
the primary in the case of X-ray ionization. The secondary 
ionization rate of H is given by 



,h _ fy-n , y(He) A He 



(18) 



and the secondary ionization rate of He is expressed sim- 
ilarly, where (0 H ) is the number of secondary ionization 
per primary electron averaged over the X-ray spectrum, for 
which we use the expression by Wolfire et al. (1995). Then, 
the terms on the right hand side mean the secondary ioniza- 
tion rate by primary electrons due to H and He ionization, 
respectively. The total ionization rate by X-rays is given by 
the sum of the primary and secondary rates: 



Cx = Cx, P + Cx,s (t = H, He). 



(19) 



In considering the secondary ionization, the H and He ion- 
ization rates are about the same mag nitudes, Cf ^ Cf e - We 
add the ionization rate by X-rays Cx (Cx) to the photoion- 
ization rate faji (£52, respectively). 



3 RESULTS 

We here present the results of our calculation and describe 
the physical processes determining the thermal evolution of 
the clouds. 



3.1 Effects of cosmic rays 

In Fig. [2] we show the temperature evolution of collapsing 
clouds irradiated by FUV fields with a diluted blackbody 
spectrum T» = 10 4 K, along with CRs with energy density 
(a) Uis = 10~ 3 , (b) 1, (c) 10, and (d) 100. The lines in each 
panel are for cases with different values of FUV strength J21 , 
which are indicated by numbers in the panel. For the cases 
with Uis = 10" 3 , 1, and 100, we also plot in Figs|3]H and 
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Figure 2. Effect of cosmic rays on the temperature evolution of primordial-gas clouds under FUV irradiation with T» = 10 4 K. 
Four panels correspond to the evolutionary tracks with four different cosmic-ray energy densities: C15 = 10 — 3 , 1, 10 and 100, where 
C15 = £/cr/10 — 15 erg cm -3 . The curves in each panel are those for different FUV fluxes, whose intensities J21 are indicated by numbers. 
The diagonal dotted lines show those at the constant Jeans mass, whose values are indicated by numbers in the Figure. 



[5] respectively, the cooling and heating rates per unit mass 
by individual processes, i.e., heating rates by compression 
(thick), by CR (dashed) and by chemical reactions (short 
dash-dotted) , and cooling rates by the H2 (thick dotted) , by 
HD (dash dotted) and by H atom (dotted). Since the major 
cooling and heating processes for [/15 = 10 are very similar 
to those for [/15 = 1, we do not repeat presenting them. 

First we see the case of C/15 = 10 -3 , where the CR is too 
weak to affect any aspect in the evolution (Figs [2] a and [3}. 
In the no-FUV case (J21 = 0), the cloud collapses along the 
standard Eb-cooling track (e.g., Palla et al. 1983). Following 
the initial adiabatic rise up to ~ 1000 K, the temperature 
starts decreasing owing to the H2 cooling, which is produced 
through the H~ channel (reactions 3 and 4), until the critical 
density ~ 10 4 cm -3 , where the H2 level-populations reach 
the local thermodynamic equilibrium (LTE) and the cool- 
ing rate per unit mass saturates. The temperature thereafter 
increases again gradually by the compressional heating to- 
wards higher density (Fig. [3] a). Addition of a FUV field 
affects the evolution in the following way. As seen in the 
cases with J21 = 0.1 and 1, with increasing radiation inten- 
sity, the initial adiabatic phase continues until higher den- 



sity and temperature, where the H2 column density becomes 
high enough for efficient self-shielding against photodissoci- 
ation. Once the H2 cooling becomes effective, the tempera- 
ture decreases and gradually converges to the track in the 
no-FUV case. For J21 = 10, before the onset of H2 cooling, 
the temperature reaches ~ 8000 K and stays almost con- 
stant around 10 2 cm -3 owing to the H Ly a cooling (dotted 
line in Fig. [2] a). This short isothermal phase is followed by a 
rapid temperature drop by the H2 cooling (see also Fig.[3]b). 
As in the cases with lower FUV fields, the temperature con- 
verges towards the no-FUV track. When a FUV field exceeds 
a critical value J21 J cr it — 20, the temperature evolution 
is qualitatively altered (top solid line in Fig. [2] a; see also 
Fig. |3]c). The isothermal evolution at ~ 8000 K continues 
until very high density (~ 10 1B cm -3 ) and the H2 cooling 
never becomes important. This bifurcation of thermal evo- 
lution originates from the fact that if the H2 formation is 
prevented until the critical density for LTE, sufficient H2 
never forms because, at higher density, (i) the amount of 
H2 needed for cooling increases (ii) collisional dissociation 
from the excited ro-vibrational levels of H2 becomes effec- 
tive, which reduces the H2 fraction (Omukai 2001). 
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Figure 3. Cooling and heating rates per unit mass as a function 
of the central number density. The cosmic ray energy density is 
Uis = 10" 3 , and the FUV intensity is (a) J 2 i = 0, (b) J 2 i = 10 
and (c) J 2 i = 20 for T, = 10 4 K. The lines show the heating 
rates by compression (thick) and chemical reactions (short dash- 
dotted), and cooling rates by H2 (thick dotted) and HD (dash- 
dotted) molecules, H atoms (dotted). 
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Figure 4. The same as Fig. \3\ but for the cosmic ray energy 
density U15 = 1, and the FUV intensity (a) J21 = 0, (b) J21 = 10 
and (c) J21 = 30. 



Next, we see how the thermal evolution described above 
changes with increasing CR flux. The CR effects are twofold, 
i.e., heating and ionization. The stronger the CR flux is, the 
more rapid is the initial temperature-increase owing to the 
CR heating as seen in Fig. [2(b) -(d). In Fig.[5]for U15 = 100, 
very high CR-heating rate can be seen at the lowest density. 
Enhanced ionization degree facilitates the H2 formation via 
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Figure 5. The same as Fig. [3] but for the cosmic ray energy den- 
sity U15 = 100, and the FUV intensity (a) J 2 i = 0, (b) J 2l = 10 
and (c) J21 = 100. In addition to cooling and heating rates shown 
in Fig. [4] heating rate by cosmic rays is also shown (dashed) . 



H~ channel. To quench H2 cooling totally, higher FUV flux 
is necessary, namely, the value of the critical FUV flux J cr i t 
becomes larger. For example, as seen in Fig. [2](d), under a 
very strong CR field of Uis = 100, even such strong FUV 
flux as J21 = 50 is not enough to quench the H2 cooling. In 
Fig- E] (a), we plot the critical FUV flux J cr it as a function 



of the CR intensity. In the case of T* = 10 4 K, from its 
low-CR limit of J cr i t — 20, Jcrit begins to increase around 
U15 ~ 10 and continues increasing as J cr it oc Uis 1 ^ 2 in the 
high U15 limit. The dominant cooling and heating processes 
are similar in all cases of the FUV flux exceeding the critical 
value Jcrit (see Figs [3] c, g] c and[S]c). 

For low FUV flux J 2 i < 0.1— 10, the higher H 2 concen- 
tration by CR ionization results in the lower temperature 
than in the no-CR case (see Fig. [5]b— d). This low temper- 
ature environment (< 150 K) favors HD formation, and its 
cooling causes further temperature decrease to a few 10 K 
(e.g., Fig. 0b-d for J21 = 0). As seen in Fig. g] (a), HD 
becomes the main cooling agent at high density and low 
temperature (n H > 10 3 cm -3 and T < 150 K). Without ex- 
tra ionization, enough HD for cooling is not produced under 
FUV radiation with J21 > 10~ 2 for T* = 10 4 K (Yoshida, 
Omukai & Hernquist 2007; Wolcott-Green & Haiman 2011). 
As seen in Figs[5](b) and[J|](b), (c), for the case of U15 = 1, 
the FUV flux higher than J 2 i ^ 1 suppresses the HD cool- 
ing and the thermal evolution becomes the same as that for 
Uis — 10~ 3 . However, the presence of stronger CRs permits 
HD cooling despite of such a high FUV flux as J21 > 10 (see 
Fig. b). 

So far, we limit our analysis to the radiation spectra 
of T, = 10 4 K. We also studied the case of T* = 10 5 K. 
As mentioned in Section 2.4.1, the H~ photodissociation 
rate depends sensitively on the brightness temperature T*. 
Then, the concentration of H2 produced through the H~ 
channel (reactions 3 and 4) changes with the value of T*. 
This leads to the enhancement of the critical FUV flux J cr it 
for T, = 10 5 K compared to the case for T» = 10 4 K (see 
Section 4). The critical flux J cr it for T» = 10 5 K is also 
plotted in Fig. g] (a). We find that J cr it = 1.6 x 10 4 for 
T* = 10 5 K, while J cr it = 20 for T» = 10 4 K in the no- 
CR cases. For CR energy density higher than C/ig ~ 10, the 
critical FUV flux J„it increases as oc Uis 1 ^ 2 in both T» = 10 4 
and 10 5 K cases. 

3.2 Effects of X-rays 

Here, we briefly mention the cases of the clouds irradiated 
both by FUV and X-rays. The temperature evolution is 
shown in Fig. □ for (a) J x ,2i = 10" 6 and (b) 10" 2 . For 
the X-ray intensity as low as Jx,2i = 10 -6 , the effects of 
X-rays, either ionization or heating, are not important and 
the evolutionary tracks in Fig. [7] (a) are the same as those 
with negligible CRs (U15 = 10" 3 ) in Fig. [2] (a). When the 
X-ray intensity is elevated to Jx,2i = 10~ 2 , the gas is heated 
instantaneously (see Fig. [7]b). In the case of J21 = 0, the 
H2 cooling balances the X-ray heating at ~ 5000 K and 
the temperature begins decreasing thereafter. For J21 ^ 0.1, 
the temperature reaches at ~ 8000 K and remains almost 
constant for a while by the H atomic cooling until H2 is 
self-shielded against the FUV field. Similar to the case of 
CR ionization, X-ray ionization enhances the critical FUV 
flux J„it for quenching the H2 cooling. This critical FUV 
flux J C rit is plotted in Fig. [S] (b) as a function of the X- 
ray intensity Jx,2i. As in the case of CR ionization, J cr it 
remains constant below a threshold around Jx,2i ~ 10~ 3 
and increases as (Jx.21) 1 ^ 2 for higher X-ray intensity. Since 
the behaviors of temperature evolution for different FUV 
fluxes as well as the dominant cooling/heating processes are 
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Figure 6. (a) The dependence of the critical FUV flux J cr it on the CR energy density C15 (solid curves) for two spectral types (thin 
curve: T* = fO 4 K, thick curve: T* = 10 5 K). In both cases, J cr it <x t/15 1 / 2 for U15 > 10. The lines indicate the relations between 
J21 and (7i5 emitted from a star-forming galaxy (equation 136 I t with (i) the Salpeter IMF and Z = 10~ 3 Zq (Pop II; thin dotted), (ii) 
Salpeter IMF and Z = (Pop III; thick dotted), (iii) the top-heavy IMF with moB = 40 Mq (short dashed) and (iv) the top-heavy 
IMF with 120 Mq (long dashed), respectively. The intersection of the curves and lines for each T* case corresponds to the actual CR 
energy density and critical FUV flux J cr it in a halo affected by the star-forming galaxy, (b) the same plot for the X-ray ionization. For 
both values of T», J cr i t oc 

(Jx,2i) 1/2 for J x ,2i > 10" 3 . Note that the fraction of halos with > 10 8 M at z = 10 irradiated by FUV 
radiation J21 exceeding > 10 5 is negligibly low < 10 — 28 . 
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Figure 7. Effect of X-ray ionization on the temperature evolution of metal-free clouds under FUV irradiation with T„ = 10 4 K. Two 
panels show the evolutionary tracks with two different X-ray intensities Jx,2l = 10 — 6 and 10 — 2 . The curves in each panel are those for 
different FUV fluxes, whose intensities J21 are indicated by numbers. The diagonal dotted lines show those at the constant Jeans mass, 
whose values are indicated by numbers in the Figure. 



very similar to the cases with CRs, we do not repeat their 
description here. 



4 DEPENDENCE OF THE CRITICAL FUV 
FLUX ON CR/X-RAY INTENSITY 

In Section 3, we see that strong CR and X-ray fluxes lead to 
the enhancement of the critical FUV flux as J cr i t oc 
and oc (Jx,2i) 1,/2 , respectively (see Fig. [6J. Also, Jcrit de- 
creases with lowering FUV temperature T». In this section, 
we discuss the reason for these dependences. 

In the cloud under a high enough FUV field, the tem- 



perature reaches ~ 8000 K and remains constant by the 
H Lya cooling. If sufficient H2 to cool forms during this 
isothermal collapse, the temperature drops by its cooling 
and the isothermal evolution is terminated. Otherwise, the 
isothermal evolution continues until very high density of 
~ 10 16 cm -3 , where the cloud becomes optically thick to the 
H~ bound-free absorption (Omukai 2001). The H2 concen- 
tration needed for its cooling to overcome the compressional 
heating is 

(3/2)fc B T 

y C ool(n2j — — — , (20) 

i-H 2 tff 
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where £h 2 = Ah 2 /^(H2) (erg s _1 ) is the cooling rate per 
an H2 molecule. Since most hydrogen is in the atomic form 
during the isothermal collapse, we assume n(H) = «h in this 
Section. In Fig. [8] we show t/ C ooi(H2) (thick solid line) as a 
function of density for the isothermal collapse at 8000 K. 
Since £h 2 increases linearly with number density until the 
critical density for LTE n H , cr ~ 3x 10 3 cm" 3 (at 8000 K) and 
saturates for higher density, j/ C ooi(H2) decreases as oc n^ 1 ^ 2 
until riH.cr, reaches the minimum ~ 1.4 x 10~ 7 there, and 
increases as oc nj/ 2 for higher density. This j/ CO oi(H2) is to 
be compared with the actual amount of H2 present ^(Hy. 
Before the temperature reaches ~ 8000 K, the H2 concentra- 
tion attains the equilibrium value set by the balance between 
the formation and destruction reactions. The effective rate 
coefficient of H2 formation via H~ channel taking account 
of the H~ photodissociation is 

k^nn 



where 



(21) 



(22) 



by using lg = 2 x 10 -7 (1 x 10 -11 ) for T» = 10 4 K (10 5 
K, respectively). This H2 is destructed either by photodis- 
sociation (reaction 18), which is dominant at the low den- 
sity where the FUV radiation is not shielded, or collisionally 
(reaction 7) at higher density. Equating the larger of those 
dissociation rates with the formation rate, we obtain the H2 
concentration 



J/(H 2 ) = min(nHfci8 >kf )teormV(e~) 
= min(j/pd(H 2 ),2/cd(H 2 )). 
The photodissociation rate coefficient is written as 



(23) 



12 (1.3 x 10- 12 ) 



(24) 

for T, = 10 4 K (10 5 



where /3is = 4.2 x 10" 
K, respectively). In the above, we assume the cloud is trans- 
parent to the photodissociating radiation since y p d(H2) im- 
mediately exceeds y c d(H2) once the self-shielding becomes 
important. The values of 2/(H2) by equation ()23|) are plot- 
ted in Fig. [8] for some combinations of T, and J21 for the 
ionization degree y(e~) — 4 x 10 -5 , corresponding the no- 
CR/X-ray case (see later). The almost vertical portion of 
2/(Ha) on the low-density side is limited by the photodissoci- 
ation (i/pd(H2)), while the gradually changing part at higher 
density is set by the collisional dissociation (y c d(H-2))- In the 
cases with J21 = 0.1 and 1 (10 3 and 10 4 ) for T» = 10 4 K (10 s 
K, respectively), y(H.2) exceeds y C ooi(H2) at some density. At 
this moment, the actual temperature falls from ~ 8000 K 
and the isothermal collapse terminates. On the other hand, 
in the case with J21 = 20 (2 x 10 4 ) for T, = 10 4 K (10 5 
K, respectively), y(H2) always remains below j/ C ooi(H2) and 
thus the isothermal collapse continues. 

For density higher than the critical density nn,cr ~ 
10 3 cm -3 , the H2 concentration decreases owing to effec- 
tive collisional dissociation from excited ro- vibrational levels 
and has no chance to reach j/ C ooi(H2) anymore. To initiate 
efficient H2 cooling, H2 concentration must exceed j/ C ooi(H2) 
before nu,ci- Thus, we can find the critical flux J cr i t by the 
condition whether J/(H2) at riH.cr ~ 10 3 cm -3 is higher than 
ycooi(H 2 ). 

In the case of irradiation with hard FUV spectra (as 
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Figure 8. The H2 concentration y(H 2 ) estimated by equation 
(123 I t (thin lines) vs. that needed for its cooling to exceed the 
compressional heating y coo i(H2) (thick line) in the isothermally 
collapsing clouds at 8000 K with the ionization degree y(e~ ) = 
4 X 10 — 5 . The H 2 concentrations are shown for the FUV fields 
with brightness temperatures T* = 10 4 K (thin solid lines) and 
10 5 K (thin dashed lines). For each value of T*, three cases with 
different FUV intensities J21 , which are indicated by numbers in 
the Figure, are presented. The dotted line shows j/ c d(H2) without 
H — photodissociation. 



T» = 10 J K), photodissociation limits the H2 concentra- 
tion and H~ photodissociation is relatively irrelevant (i.e., 
fc form ~ k 3 ). From the condition j/pd(H 2 ) = y C ooi(H 2 ) at 
?iH,cr, we obtain 



J crit 



1 fe 3 W H ,crj/(e 

J/cool (H2) 



(25) 



On the other hand, for softer FUV spectra, (as T» = 10 4 
K), the H 2 formation rate is significantly reduced by H~ 
photodissociation: fcf orm = k^kinn/kiQ. The dotted line in 
Fig. [8] shows j/cd(H2) without H~ photodissociation, which 
exceeds j/ C ooi(H 2 ) before reaching riH.cr- The actual value 
of j/cd(H 2 ) is, however, suppressed by IT" photodissociation 
(see thin solid lines in Fig. [8} and thus the H 2 concentration 
is limited by collisional dissociation. Therefore, we obtain 
from the condition y c d(H 2 ) = y C ooi(H 2 ) at n H ,c r : 



J c] 



ki fc 3 n H , cr ^(e ) 
Pi9k 7 j/ coo i(H 2 ) 



(26) 



In general, the critical FUV flux is given by the larger of the 
above two values (equations 1251 and \2Q \ . 

Note that J cr i t increases with the ionization degree 
y(e~), which can be evaluated as follows. Without external 
ionization, the ionization degree is governed by 



— — — = -~k 2 n H y{e ) . 



(27) 



Integrating this equation for a collapsing cloud (i.e., equa- 
tion [[} , 



V(e ) = 



3/0 (e ) 



1 + 2fc22/o(e-)n H iff ' 



(28) 



where yo(e ) is the initial ionization degree. For density > 
50 cm -3 , this renders y(e~) = (2k2nntg)~ 1 . Equation 
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gives y(e ) ~ 4 x 10 5 at nn = 
at 10 3 cm -3 , equations (|25|) and 



Using the values 
lead to the critical FUV 
flux Jcrit ~ 10 (10 4 ) for T* = 10* K (10 s K, respectively), 
which are indeed similar to our numerical results for the 
no-CR/X-ray cases (Fig. [6]). 

In the presence of strong CRs, the ionization degree is 
set by the balance between recombination and CR ioniza- 
tion: 



CCRA 1 / 2 -1/2 



(29) 



Comparing equations (|28l) and (|29[) . we see that the CR 
ionization controls the ionization degree if its rate is higher 
than 



CcR.crit — 



3.8 x 10" 19 s"\ 



(30) 



where the last equation is evaluated at T = 8000 K. This 
value is equivalent to the CR energy density (TcR.crit — 5.5 x 
10~ 15 erg cm -3 from equation (fl2|) . For Ucr > E^CR.crit, the 
ionization degree, as well as the critical FUV flux, increases 
as oc Ucr 1 ^ 2 - This relation is reproduces well the behavior 
of Jcrit as a function U15 shown in Fig. [!( a). 

Note that our critical CR density is similar to another 
threshold value where the HD cooling becomes important in 
the absence of FUV radiation, Ucr = 2 x 10 -15 erg cm -3 
(Stacy & Bromm 2007). This is not a coincidence: both crit- 
ical values are related to enhancement of the ionization de- 
gree and thus H2 concentration from the case without CR 
ionization. 

Next we consider the case of the X-ray ionization. Just 
the same as in the CR ionization, X-rays enhance the ioniza- 
tion degree as y(e~) oc (Jx,2i) 1,/2 and thus J cr it oc (Jx^i) 1 ^ 2 
for Jx.21 > 10 -3 , whose threshold is determined by the bal- 
ance between recombination and X-ray ionization. This de- 
scribes well the results shown in Fig. [6](b). 

Our value of J cr it = 20 for T» = 10 4 K is smaller than 
./crit = 39 which Shang et al. (2010) estimated with one- 
zone model. This difference comes from the liVformation 
rate at high temperature: the rate coefficient we adopt (Galli 
& Palla 1998) is smaller than that Shang et al. (2010) used 
(Abel et al. 1997). On the other hand, our value of J cr i t = 
1.6 x 10 4 for T* = 10 5 K is slightly larger than Shang et al. 
(2010) because we use the new shielding factor by Wolcott- 
Green & Haiman (2011). 



5 CR AND X-RAY FROM STAR-FORMING 
GALAXIES 

So far, we have regarded the intensities of FUV, CRs and 
X-rays as free parameters in our calculation. All those radi- 
ation fields are, however, closely linked to the star-formation 
activity in galaxies and thus their intensities, i.e., J21, U15 
and Jx,2i, are all proportional to the star formation rate 
(SFR). We here consider the star-forming galaxies with two 
types of initial mass function (IMF): the Salpeter IMF and 
the top-heavy one. For the Salpeter IMF, we take the mass 
range of 1 — 100 Mq and consider two cases of the stellar 
metallicity Z = 10 Zq and 0, which corresponds to Pop 
II and III star clusters, respectively. For Pop III galaxies, we 
also consider the two cases of the top-heavy IMF where all 



Table 2. The IMF models of star- forming galaxies 



model 


IMF 


mass (Mq) 


metallicity (Zq) 


r, (K) 


(i) 


Salpeter 


1 - 100 


10~ 3 


10 4 


(ii) 


Salpeter 


1 - 100 





10 s 


(iii) 


top-heavy 


40 





10 5 


(iv) 


top-heavy 


120 





10 5 



stars are m B = 40 Mq or 120 Mq. The IMF models we 
study are summarized in Table 2. 

For FUV flux, we use the Lyman- Werner photon emis- 
sivity from a star-forming galaxy for the IMF models calcu- 
lated by Schaerer (2002; 2003): 



J21 



6.5 x 10 2 

1.6 x 10 3 
1.4 x 10 4 



d 



SFR 



-), (31) 



,10 kpc/ V20M Q yr 

1.7 x 10" 

where d is the distance from the source, and the numbers 
in the bracket correspond to the IMF models in the same 
order as in Table 2. 

We assume that supernova remnants (SNRs) are the 
sources of CRs. The total CR energy from a SNR is written 
as (Stacy & Bromm 2007) 

-Esn 



E C r = 10 



m(PCR\(_ 

\ 0.1 / V: 



erg, 



(32) 



, 10 52 erg/ 

where pgr is the fraction of SN explosion energy _Bsn going 
into the CR energy. Using the supernova rate in the source 
galaxy 

foB SFR 



iVsN 



moB 



(33) 



where /ob, mob is the mass fraction and average mass of 
massive (^8 Mq) stars, the CR energy density of Ucr is 

Ecr 



Ucr = iVsN- 



(34) 



47rd 2 (v)' 

where (v) is the average CR velocity. For the fiducial values 
of Pcr = 0.1 and Esn = 10 52 erg, we obtain 



Uis = 4 



42 
42 
44 
15 



, 10 kpc 

From equations (|31[) and (I35|) 



d 



H 



SFR 



J21 = < 



15 

38 

3.1 x 10 2 
1.4 x 10 3 



20 Mq yr- 



> Ui 



(35) 



(36) 



X-rays are mainly emitted by high-mass X-ray binaries. 
Observationally, the X-ray luminosities in 2 — 10 keV of local 
star-forming galaxies are related with their star formation 
rates as (e.g., Glover & Brand 2003) 

, ( SFR \ _ a 



Lx = 1.2 x 10 4 



namely, 

Jx,2i = 4.5 x 10" 



v20 Mq yr- 



d 



10 kpc 



SFR 



20 Mq yr- 



(38) 
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Although this is an empirical relation for local galaxies, ob- 
servations of Lyman-break galaxies provide supports for this 
to be valid as far as z ~ 4 (Glover & Brand 2003 and ref- 
erences therein). Here, we thus extend the relation (|38[) to 
high-redshift (say, z ~ 10) universe. From equations (|31[) 
and (|38[l . we obtain 



J21 = < 



1.4 x 10 J 

3.5 x 10 5 
3.0 x 10 6 
3.7 x 10 6 



> Jx.: 



(39) 



In Fig. |6] (a), we plot the relations J36j for the IMF 
cases (i)-(iv) in Table 2. The FUV intensity and CR energy 
density from galaxies are expected to fall on this relation. 
In the case of Pop II radiation sources (i.e., type i), which 
are characterized by T* ~ 10 4 K, the relation (|36|1 intersects 
with the Jcrit curve at (J21, U15) = (20, 2), where the J cr it 
does not deviate yet from its no-CR value. This means that 
the CR ionization does not modify the value of J cl -i t signifi- 
cantly. On the other hand, for Pop III star clusters (T» ~ 10 
K), the critical FUV flux at the intersection points are > 20 
(type ii), 10 (iii), and 2 times (iv) as large as J cr i t in the no- 
CR case, respectively. Namely, the CR ionization drastically 
changes the critical FUV flux except for the case with the 
top-heavy IMF with m OB = 120 Mq. 

In Fig. [6](b), we also plot the relations t|39|l between J21 
and Jx,2i • As in the case of CRs, the X-ray ionization effects 
become important for higher T„, in particular, T* > 10 4 K. 



6 CONCLUSION AND DISCUSSION 

We have calculated the thermal evolution of primordial 
clouds under strong far-ultraviolet (FUV) fields, along with 
cosmic rays (CRs)/X-rays. In the cloud under a FUV field 
exceeding a threshold, the H2 cooling is suppressed at any 
density. Such a cloud collapses almost isothermally at 8000 
K by hydrogen atomic cooling. According to numerical sim- 
ulations (Bromm & Loeb 2003), it avoids fragmentation and 
collapses directly to a supermassive star (SMS). 

Without external ionization by CRs/X-rays, the crit- 
ical FUV flux is J cr it ~ 1.6 x 10 4 (in units of 
10 -21 erg s _1 cm -2 sr _1 Hz -1 ) for a diluted black body 
spectrum with brightness temperature T» = 10 5 K, while 
it is as small as J21 — 20 for T, = 10 4 K. This dependence 
on T* comes from the higher H~ photodissociation (> 0.755 
eV) rate for lower T* at the same J21 normalized at 13.6 eV. 
Since H2 is produced through the H~ channel (reactions 3 
and 4 in Table 1), the larger H~ photodissociation reduces 
the amount of H2 and thus the value of J cr it . 

We have studied how this critical FUV flux changes 
with simultaneous irradiation of either CRs or X-rays. In 
the case of CR irradiation, the critical FUV flux J cr i t 
begins to increase with the CR energy density Ucn ^ 



10~ 14 erg cm -3 and depends as J cr i t oc Ucr 1 ^ 2 asymp- 
totically. Similarly, for X-ray irradiation, J cr it oc Jx 1//2 for 
Jx > 10 -24 erg s _1 cm -2 sr _1 Hz -1 at 1 keV. In both cases, 
the increase of the critical FUV flux is caused by the en- 
hanced ionization by CRs/X-rays, which promotes the H2 
formation and cooling. 

FUV intensity and CRs/X-rays from star- forming 
galaxies are expected to correlate each other since they all 



trace the massive-star forming activity. Using the expected 
relations between FUV intensity and CRs/X-rays, we have 
found that if the initial mass function (IMF) of the radia- 
tion source is Salpeter-like and the brightness temperature 
of the FUV radiation is rather high T„ ~ 10 5 K, the criti- 
cal FUV intensity increases significantly J cr i t ~ 10 6 owing 
to ionization by CRs and X-rays. Even with the top-heavy 
IMF, the critical FUV intensity increases to J cr it ~ 10 s un- 
less the stellar mass is > 100 Mq. Since the fraction of halos 
exposed to FUV flux exceeding J cr j t decreases exponentially 
with J cr it for > 2 x 10 3 (Dijkstra et al. 2008), there is little 
possibility (< 10 -28 ) for such intense FUV field realized in 
any halos. We conclude that if the radiation source is com- 
posed of Pop III stars with brightness temperature ~ 10 5 K 
and contains sources of CRs/X-rays, its IMF must be very 
top-heavy > 100 Mq to enable SMS formation in nearby 
halos. Other possible radiation sources enabling SMS forma- 
tion include, e.g., sources with low brightness temperature, 
T„ ~ 10 4 K, such as Pop II /I star clusters, or those with 
high T* but without strong CRs and X-rays. The latter pos- 
sibility is realized if a Pop III star cluster is so young that 
it harbors neither supernova remnants nor high-mass X-ray 
binaries. 

We have also found that the extragalactic ionization ef- 
fects (CRs/X-rays) are able to promote the HD formation 
and cooling even if the clouds are irradiated with FUV ra- 
diation. Without external ionization, Yoshida et al. (2007) 
and Wolcoott- Green & Haiman (2011) showed that the HD 
formation is suppressed by FUV radiation with J21 > 10 -2 
for T„ = 10 4 K, and that the number of metal-free stars 
(> 10 Mq) formed in Hn regions by HD cooling (so-called 
Population III. 2 stars) is reduced. However, we find the pos- 
sibility that such stars can form under FUV radiation if 
strong CRs/X-rays ionize the clouds. These mean that the 
star formation in Hn regions by HD cooling is not suppressed 
so strongly as previous thought (see also Wolcott- Green & 
Haiman 2011). 

Finally, we discuss limitations of our analysis. We have 
treated the cloud evolution at the center by the one-zone 
model, where two assumptions have been made. First one is 
that the density increases approximately at the free-fall rate 
(equation [TJ. This is a good approximation for the isother- 
mal spherically symmetric collapse. In reality, however, in 
addition to the thermal pressure, the turbulent support and 
compression, as well as the centrifugal support, which we 
have not taken into account, can have importance in dy- 
namics, which controls the collapse rate. For example, in 
relatively massive halos, which we have considered here, the 
primordial gas is known to be in a state of turbulence (Wise 
& Abel 2007; Greif et al. 2008). Such turbulence might affect 
the thermal evolution, as well as the exact value of the FUV 
critical intensity for suppressing H2 cooling. If the angular 
momentum is present, the centrifugal force becomes more 
important with contraction relative to the gravity, and even- 
tually halts the collapse. We here, however, only considered a 
rather low-density regime, and thus the angular momentum 
effect is not so dynamically important in our case. Another 
assumption in our model is that the external radiation is 
attenuated with the column density estimated by the cen- 
tral density and the size of the core, which is given by the 
Jeans length (equation [7]). Although the column density of 
the core is given by this value for the Larson-Penston solu- 
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tion, the contribution from the envelope needs to be added. 
But we can see the envelope contribution is not significant 
by the following consideration. The steep decline of density 
with radius as p oc r~ 2 limits the column density of hy- 
drogen nuclei A^h in the envelope at most that of the core. 
Due to the photodissociation and low formation rate, the 
H2 abundance is remarkably lower in the envelope than in 
the core. Therefore, the envelope contribution to the col- 
umn density of H2, which is relevant for evaluating the crit- 
ical FUV intensity, is small and the H2 column density is 
well approximated by that of the core. Another concern is 
about the geometry of the cloud. In the low-density regime 
we considered, the cloud shape can strongly deviate from 
the sphere. Even though the clouds are in such shapes as 
sheet-like or filamentary, however, the length scale of the 
shortest axis, which is most relevant for the shielding ef- 
fects, is still roughly given by the Jeans length, and thus 
the our assumption of the Jeans length shielding remains 
valid. In fact, Shang et al. (2010), who studied the evolution 
of the primordial clouds under strong FUV irradiation and 
compared the results by the one-zone model and those by 
the cosmological three-dimensional simulation, found very 
good agreements as for the thermal evolution at the center 
of the clouds, as well as for the critical FUV intensity J cr it 
needed to quench the H2 formation/cooling. In addition, we 
assume that the clouds are isotropically irradiated by FUV, 
as well as CRs and X-rays. However, such intense FUV field 
as exceeding the critical value J cr i t tends to be dominated 
by a nearby single large source, rather than by a collective 
effect of a large number of sources (Dijkstra et al. 2008). We 
suspect that, in the case of such anisotropic radiation field, 
more intense sources are necessary to induce SMS forma- 
tion as the radiation comes from only limited solid angles. 
To confirm it, however, detailed modeling is required, which 
is beyond the scope of this paper. 
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